Fluctuation-induced drift in a gravitationally tilted optical lattice 
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Experimental and theoretical studies are made of Brownian particles trapped in a periodic po- 
tential, which is very slightly tilted due to gravity. In the presence of fluctuations, these will trigger 
a measurable average drift along the direction of the tilt. The magnitude of the drift varies with 
the ratio between the bias force and the trapping potential. This can be closely compared to a 
theoretical model system, based on a Fokker-Planck-equation formalism. We show that the level of 
control and measurement precision we have in our system, which is based on cold atoms trapped in 
a 3D dissipative optical lattice, makes the experimental setup suitable as a testbed for fundamental 
statistical physics. We simulate the system with a very simplified and general classical model, as 
well as with an elaborate semi-classical Monte-Carlo simulation. In both cases, we achieve good 
qualitative agreement with experimental data. 

PACS numbers: 05.60.-k, 05.40.Jc, 37.10.Jk 



I. INTRODUCTION 



A very general problem in physics is that of a Brow- 
nian particle moving in a periodic potential; a seminal 
treatment, using Fokkcr-Planck formalism, is given by 
Risken [l[ . Of particular interest is the 'tilted washboard 
potential', where the Brownian particle is also subjected 
to a constant force, which can actually be used to model 
a wide variety of physical systems (sec, e.g., [2j, |3[ and 
references therein). Recently, there has been a signifi- 
cantly increased interest in the dynamics of small sys- 
tems, where fluctuations and noise play a dominating 
role and where a classical thermodynamic equilibrium 
does not occur. There have been theoretical discover- 
ies (e.g., |^-i8J) providing understanding of fluctuations 
and non-equilibrium situations, as well as experimental 
breakthroughs (e.g., |9|, [l(|). Closely related to this are 
systems where noise, or fluctuations, is the source for 
directed drift, so-called Brownian motors (see, e.g., (Ill — 
Il3j|). or where the noise opens up a possibility for drift 
in a biased system, where this bias would otherwise not 
have been enough to overcome potential barriers and/or 
friction. 

In this work, we trap and hold cold atoms for several 
seconds in a three-dimensional, dissipative optical lat- 
tice [lj, [l5j . The thermal energy of the atoms is of the 
order of a tenth of the depth of the periodic potential, 
and the tilt of the 'washboard potential' in the vertical 
direction due to gravity is approximately three orders of 
magnitude lower than the potential depth (on the range 



of the period of the potential) [16J : thus the potential 
should support the atoms from gravity with a very good 
margin. However, these dissipative optical lattices put 
the atoms in a regime where fluctuations play a domi- 
nating role, and where dissipation is also present. These 
fluctuations will trigger a discernible drift, even with such 
a small bias, and threshold effects will be present if pa- 
rameters, such as external force, potential depth, fluc- 
tuation amplitude, or damping, are varied. This makes 
this system a suitable experimental testbed for funda- 
mental studies of fluctuation phenomena. In addition, 
the setup used here is very close to the double optical 
lattice arrangement that has been used to create a Brow- 
nian motor J17H19| . The present work is therefore also 
of interest for an understanding of the role of gravity in 
that context. 



II. THE PROBLEM IN A FOKKER-PLANCK 
EQUATION CONTEXT 

A classical particle in the above predicament, with 
(vertical) position coordinate x, will follow the Langevin 
equation 

x = —- r v(x)-~rx + - + t(t). (l) 

m ax m 

Here, m is the mass, 7 is a uniform damping constant, F 
is a uniform external force, and £ is a Langevin stochastic 
force PJ. The periodic potential is 



V(x) = V sm(2irx/L) . 



(2) 
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where L is the spatial period of the potential. 

The characteristics of such a system will be determined 
by the relative strengths of the terms in Eq. (JlJ, and 
in particular the magnitude of the friction is important. 



The case that compares most closely with our system is 
one where the friction is relatively small, or in another 
terminology, where the system is not overdamped. This 
means that the particle can be either in a 'locked state', 
where it oscillates around a minimum in one potential 
well, or in a 'running state', where it travels from well to 
well; and it will undergo transitions between these states. 
The mobility of an ensemble of particles is defined as 
/x = (x)/F and, for the frictionless case, the locked and 
running states would correspond to 7/1 = and 7// = 
1 respectively. The solution to this general problem is 
outlined in p| . In the case where the friction is small but 
still significant, and the noise term is of the same order as 
the potential depth (in appropriately rescaled units), the 
transition from locked to running, with increased force 
(or decreased potential depth), is much less sharp, and 
the mobility never quite becomes zero, as long as there 
is some noise. For the case where the noise term and the 
friction are not spatially dependent, there exist analytical 
solutions to this general problem [l| . 



of the mechanisms for the continued route to equilib- 
rium, for an atom localized in a well, are not precisely 
known. However, with strong support from experimen- 
tal and theoretical investigations (c./., [22], |24H29I |). we 
here assume the following. An atom trapped in a well 
experiences no direct damping. However, its probabil- 
ity of acquiring energy from light scattering, and to get 
unlocked, is higher the more excited it is in the well. 
When it gets unlocked, it will be again be exposed to 
laser cooling, it will loose its kinetic energy, and it will 
be trapped again in some bound state. As this goes on, 
there will be a gradual accumulation towards lower lying 
and more deeply trapped states, from which the escape 
probability is low, and eventually an equilibrium will be 
reached. Furthermore, the deeper the potentials are, the 
larger the portion of atoms that are trapped, but even 
for very shallow potentials the majority of the atoms are 
trapped. Or, correspondingly, one atom spends most of 
its time being trapped, interrupted by short periods of 
inter- well flight [24|, [23, U&, bM| > where it can travel over 
several wells. 



III. A TILTED DISSIPATIVE OPTICAL 
LATTICE 



B. The damping term in the current work 



Dissipative optical lattices arise from atom interaction 
with a periodic light shift potential, created by a num- 
ber of laser beams, tuned below and relatively close to 
an atomic, dipole-allowed transition [14|,[l5|. In our case, 
the detuning, A, is typically of the order of 10-40 natural 
linewidths, T, of the transition in question. The proxim- 
ity to a resonance means that incoherent light-scattering 
will be important for the dynamics of the atoms. There 
will be diffusion effects associated with photon recoils and 
with instantaneous changes to the light shift potential, 
due to optical pumping [20J; these 'heating' effects cor- 
respond to the noise term, £, in Eq. (JTJ). Moreover, with 
a proper configuration of the laser beams, laser cooling 
will be present (c./., [lj, Qji [20, [2l| ) ; corresponding to 
the damping, 7. 



A. Laser cooling 

The seminal treatment of laser cooling [2(J, [22| is really 
only relevant for atoms that move around in the lattice, 
and the approximate approach is there taken that the 
friction constant, 7, is a constant and that a spatial av- 
erage can be used. This allows for a reasonably straight- 
forward treatment based on a Fokkcr-Planck equation 
approach [22j. However, in actual experiments with dis- 
sipative optical lattices, this cooling mechanism may be 
relevant for the initial damping of the thermal energy, 
and for the first phases of the route to equilibrium, but 
the atoms will quickly loose enough thermal energy in or- 
der to be trapped in the potential wells of the lattice; and 
at equilibrium, they indeed typically get localized close 
to the bottom of the potentials (c./., [23|). The details 



In the current work, the relevance of trapping strongly 
affects how the damping is to be treated. We make 
the working hypothesis that when an atom is trapped 
('locked state'), its motion is undamped. If, and when, 
it becomes untrapped ('running state'), an effective fric- 
tion, 7lc (with LC standing for "laser cooling"), turns 
on, which we assume can be reasonably well approxi- 
mated by the spatial average used in [201 ] . i.e., an un- 
trapped atom is subjected to dissipation of its momen- 
tum as in the traditional picture of laser cooling. 

The acceleration due to the tiny bias force, F = mg 
(where g is the gravitational acceleration), is so small 
that it will not significantly affect the velocity of the 
atoms during a single period of inter- well flight. The 
damping, — Tlc^j will occur only due to the velocity the 
atoms acquire from the light scattering, i.e., the Langevin 
force, £, and a free atom will shortly be trapped again. 
Thus, the dynamics of the atom will be of a 'stop-and-go' 
nature. The effect of gravity will only be a very small av- 
erage downward drift of the center-of-mass of the sample, 
partly because a trapped atom has a slightly higher prob- 
ability to escape downwards than upwards, and partly 
because an untrapped atom will travel slightly longer 
distances when going downhill than uphill. Thus, (x) 
is assumed to change (downwards) linearly with time, 
since any memory of the gravitational acceleration will 
be erased when the atom is recaptured (see also |29|V 



IV. EXPERIMENT 

For the vast majority of experiments done with dissipa- 
tive optical lattice, the holding time in the optical lattice 



has been rather short, i.e., 10-100 ms. This is partially 
because it is difficult to achieve longer lifetimes, and more 
importantly because when the laser cooling dynamics is 
studied, longer time scales have not been believed to be 
important. The basic idea behind our experiment is to 
hold the atoms for much longer times, approaching 10 s, 
and study how the mean position of the sample evolves. 
This can be done by direct imaging of the atoms in situ. 
However, as a more precise diagnostic, we release the 
atoms and measure their arrival at a laser probe located 
at a distance, I = 5 cm, below the sample ('time-of-flight 
detection' |lj|). We do this for a range of different poten- 
tial depths, Vq, providing us with data for the mobility 
as a function of F/Vq. The potential depth is varied by 
adjusting the irradiances and the detunings of the opti- 
cal lattice laser beams [lj] . From the time-of-flight data, 
we can also analyze the velocity distribution in more de- 
tail, in order to approximately quantify how much time 
an atom spends in the locked and the running states on 
average [31| , 



A. Experimental set-up 

The experimental set-u p h as been described in more 
detail elsewhere [18|, Qjl, bL0, [32|- In short, we trap 
and cool cesium atoms with standard laser cooling 
techniques |21| . The initial cold cloud typically con- 
sists of about 10 8 atoms with a temperature of around 
5 /xK, which corresponds to about 25-E rec , where E Tec = 
p^ ec /2m, is the kinetic energy associated with the recoil 
of absorption or emission of a single resonant infrared 
photon (of momentum p rec ). The atoms are then trans- 
ferred to a three dimensional optical lattice, which is 
constructed by four laser beams in a three dimensional 
generalization of the 'lin_l_lin configuration' |14l . Il5l |33| . 
The optical lattice is detuned by either 30 or 40 natural 
lincwidths, T (T = 2tt x 5.21 MHz), from the transi- 
tion between the hyperfine structure states F g = 4 and 
F c = 5, within the D2-line of Cs [H]. This configura- 
tion will give a phase-stable lattice with a face-centered- 
tetragonal geometry, and with a vertical symmetry axis. 
The lifetime of the atomic sample is limited by collisions 
with the residual gas in the vacuum chamber and by dif- 
fusion in the optical lattice. The maximum lifetime is of 
the order of 10 s, giving us ample time to perform the 
intended studies. 
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FIG. 1: (Color online) Position of the center-of-mass of the 
atomic cloud derived from time-of-flight detection data, as a 
function of holding time r in the optical lattice for the po- 
tential depths 45-E r oc (black), 95-E rC c (red), 125-E rcc (green), 
200-Erec (blue) and 395£ rcc (purple), for a detuning of — 40F 



the drift, evidencing the 'stop-and-go' dynamics, is ob- 
served from the data presented in Fig. [T] where a faster 
drift due to the gravitational tilt is also observed for shal- 
lower potentials. 



1. Mobility 

In order to extract the velocity of the drift (x), straight 
lines are fitted to curves like the ones in Fig. [I] Figure [2jt 
shows results for a range of data, with (i) as function of 
the potential depth for the detunings — 30r and — 40r. 
In Fig. [^b we display the same data, but now scaled as 
the mobility fiasa function of the constant force F di- 
vided by the potential depth. This allows for a more 
direct comparison with the general theoretical treatment 
in [l|]. Such a comparison must be made with care, since 
in || a spatially uniform friction is assumed, whereas 
our experimental conditions arc such that the damping 
force that acts on individual atoms depends strongly on 
position and kinetic energy. Moreover, even with a sim- 
plified model for a friction force, the coefficient of friction 
will vary with detuning [20|. Nevertheless, fluctuation- 
induced drift in the tilted potential is clearly demon- 
strated. 



B. Experimental results 

In Fig. [TJ we show measured positions of the center of 
mass of the atomic cloud as a function of holding time 
in th e op tical lattice, r, as derived from time-of-flight 
data [19|. The arrival time of an atom to the time-of- 
flight probe is a function of its initial position and initial 
velocity. From the average arrival time, the average posi- 
tion, (x) = {x)t, follows directly. The linear evolution of 



2. Running and locked states 

By analyzing the velocity distributions obtained by 
time-of-flight detection, the fraction of atoms in the run- 
ning state compared with the total amount of atoms, 
N Tun /Ntot , can be extracted [19|, (27|, [31j . Assuming that 
the momentum distribution corresponds to a (truncated) 
Gaussian core of trapped atoms, with wide wings corre- 
sponding to untrapped atoms, we can calculate approxi- 
mate numbers for N Tun /N tot and by Gaussian fits to the 
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FIG. 2: (Color online) (a) Derived drift velocity, {£), as a 
function of potential depth, Vb, for the detunings — 30F (open 
squares), and — 40F (circles), (b) the same data, but plot- 
ted as the mobility, (iasa function of (FL)/Vo- The solid 
line is from a semi-classical Monte-Carlo simulation, while 
the dashed line is from a simplified classical simulation. 
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FIG. 3: (Color online) (a) The fraction of unlocked atoms, 
as a function of potential depth, Vb, for the detunings, — 30r 
(open squares), and — 40F (circles) The solid line represents 
the semi-classical Monte-Carlo simulation, and the dashed 
line is the classical simulation, (b) The drift divided by the 
fraction of atoms is the running state versus the potential 
depths. A straight line indicates the recoil velocity. Note 
that velocity is defined to be positive upwards 







_._— 














_^J^-~ 




D 






• >^l 




• 
a 




D 




/B D 




• 








p: u 

I? 


• 






















• /.. 












;• 













SIMULATIONS 



Semi-Classical Monte-Carlo Method 



momentum distributions. The results of this for the same 
data used to observe the drift is shown in Fig. [3^. Within 
our range of detunings and irradiances, the fraction of the 
atoms that are free is typically just a few percent. For 
potential depths smaller then about 200-E roc , the fraction 
of atoms in the running state increases drastically, and 
for the shallowest potentials we use, it gets as high as 
about 25%. This behaviour has striking similarities with 
Fig. [2^, where the average velocity downwards due to 
gravity also increases drastically for the same potential 
depths. As an attempt to investigate the dependence be- 
tween the drift and the fraction of atoms in the 'running 
state' further, in Fig. [3p we plot {x)/(N Iun /N to t) versus 
the potential depth. From this it is clear that the drift, 
not being constant, depends not only on the fraction of 
untrapped atoms but also on something else. 



With a semi-classical Monte-Carlo simulation of the 
atom- laser interaction [30, |35( , we derive theoretical data 
corresponding to all experimental curves shown. In these 
simulations, the laser field and the motion of the atoms 
are treated classically, which enables tracking of the po- 
sition and momentum of each particle, while the in- 
ternal state of the atom is treated quantum mechan- 
ically, using the true degenerate level structure of the 
F g = 4 — > F c = 5 transition. In addition, the presence of 
the excited F c — 4 level is also included [35|. Diffusion 
and friction arise "naturally" from the laser-atom inter- 
action and are position and velocity dependent, as in the 
experiment. 

The simulations comprise 15000 non-interacting 
atoms, starting from an initial cloud contained in a single 
lattice well, at a temperature of 5 /iK. The drift veloc- 



ity is obtained by calculating the median position of the 
atoms at the end of the simulation and dividing by the 
duration of the simulation (<~ 50 ms). This is necessary 
because a direct calculation of (i) is affected by the pres- 
ence of a few high-velocity atoms, corresponding atoms 
that become untrapped and never get recaptured by the 
lattice. Such atoms arc believed to be also present in 
the experiment, but drift out of the lattice and never get 
detected. 



The optical lattice is here one dimensional (along the 
vertical), which means that an exact quantitative agree- 
ment cannot be expected, but all qualitative features in 
the experiment are reproduced, and so are the orders 
of magnitude of the mobility, the fractional populations, 
and the potential depths where significant features oc- 
cur. Also, the fraction of atoms in the running state is 
calculated much more precisely from total energies and 
positions of the individual atoms, than it can be deter- 
mined from the experimental time-of-flight data. The 
results, smoothed out to remove the fluctuations due to 
the small sample size, are shown together with the ex- 
perimental data in Figs. [5] and [3] 

One of the constraints of the current experimental 
setup is the use of gravity as the external force. This 
means that the variation of FL/Vq can only be achieved 
by a change in Vq. However, Vq is not independently con- 
trollable, but it depends on the irradiance of the lasers 
and their detuning with respect to the atomic transition. 
These in turn affect friction (Sisyphus cooling) and diffu- 
sion (photon scattering), and therefore also the temper- 
ature of the atoms, such that it can be seen as being a 
function of the optical lattice potential depth (see, e.g., 
011). 
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FIG. 4: (Color online) Mobility calculated from semi-classical 
Monte-Carlo simulations. +: constant force F = mg and 
varying Vb [same data as full line in Fig. HJb)]; o: varying 
F and constant Vb = 7QE lac ; X: varying F and constant 
Fo = U0E rcc . The detuning is -30T. Panels (a) and (b) 
differ only by the scaling of the abscissa. 



B. Classical Approach 



To investigate how this affects the mobility, we have 
run simulations at constant potential depth for different 
values of the external force. In Fig. Sta), the simulation 
result previously shown in Fig. [UJb) is plotted along with 
the mobility obtained varying F, for two different values 
of Vq. Over the range covered by the experiment, the 
mobility now appears constant, with a value dependent 
on the potential depth. Further increasing the value of 
the external force, Fig. Sfb), the well-known behavior of 
Brownian particles in tilted potentials emerges, with a 
transition region between locked and running states [J. 
It thus appears that the variation of the mobility as a 
function of Vb, Fig. [51 reflects the fact that it is the laser 
irradiance, and not strictly the potential depth, that is 
modified. In addition, the current experiment probes a 
region corresponding to a low-mobility locked state, in 
accord with the fact that most atoms have an energy 
below the well-to- well barrier height [23, bL!| . Note that 
the semi-classical model used here cannot reproduce the 
Dopplcr cooling [20|, [2l| that will become important as 
the velocity of the atoms increases. 



In order to check the generality of the system and rele- 
vance of an analysis in terms of a Fokker-Planck equation, 
as in Eq. (TT|), we also perform a simple, completely clas- 
sical Monte-Carlo simulation for the Langevin equation 
corresponding to Eq. (p}, for a ID system with classical 
particles in a tilted washboard potential, as in Eq. @. 
In this simulation, we let the noise term and the friction 
scale linearly with the potential depths. For simplicity 
and generalization we take a uniform friction. The re- 
sults are presented together with the experimental data 
in Figs. [2] and [3] The basic characteristics of our system 
are reproduced. 



VI. CONCLUSION 

To conclude, we have made a quantitative study of 
how random isotropic fluctuations, together with a very 
small bias force, gives rise to an average drift. When the 
fraction between bias force and trapping potential, F/Vq 



is varied, the magnitude of this drift can change. The 
system can be well described by a Fokker-Planck equa- 
tion formalism, and the experimental control, together 
with the precision in the measurements, make the sys- 
tem suitable as a general testbed for studies of funda- 
mental fluctuation phenomena. To emphasize this, we 
qualitatively reproduce our data with a simplified clas- 
sical simulation, as well as with a careful semi-classical 
Monte-Carlo simulation of the laser cooling setup. Our 
results also evidence the 'stop-and-go' nature of the dy- 
namics of the atoms, where they continuously exchange 
between being t rapp ed in p otential wells and travelling 
over many wells II E3, Hi, H3] ■ 

One of the constraints of the present experiment is the 
use of gravity as the bias force. This leaves the potential 
depth Vq as the main variable parameter, but it is only 
accessible through the laser irradiance, which also mod- 
ifies diffusion and friction. One possible solution would 



be the introduction of an additional laser beam, using 
its radiation pressure as the bias force. Moreover, it is 
possible to operate the optical lattice far-detuned from 
the atomic resonance, where it only serves as a (conserva- 
tive) potential. Extra laser fields could provide diffusion 
and friction, allowing the exploration of a broad range of 
scenarios. 
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